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BLAST  CHARACTERISTICS  OF  20  and  100  TOR  HEMISPHERICAL 
AH/FO  CHARGES,  NOL  DATA  REPORT 


Prepared  by: 

1.  D.  Sadvin 
K.  M.  Svisdek,  Jr. 


ABSTRACT:  Two  twenty-ton  and  one  100 -ion  Hemispherical  AN/FO  (toaoniua  nitrate/ 
fuel  oil)  charges  were  detonated  on  the  surface  at  the  Defence  Research  Esteblishmenl 
Suf field,  Ralston,  Alberta,  Canada.  The  testa  were  conducted  during  A>igust  1^69 
aa  a  cooperative  U.S. /Canadian  effort. 

The  major  results  were: 

1.  AH/FO  hee  been  demonstrated  to  be  a  highly  suitable  explosion  source  for 
nuclear  airblaat  simulation. 

2.  Over  the  1-200  pgi  region,  there  was  co  significant  difference  iu  the 
pressure-distance  characteristics  between  AN/FC  and  TNT. 

3-  The  impulse  characteristics  of  the  AN/FO  system  were  found  to  be  slightly 
lower  than  those  of  TNT. 

h .  Ho  self  heating  of  AN/PO  was  observed. 

3 

5.  Conventional  cube  root  scaling  applies  for  AN/ TO  over  a  10  range  in 
explosive  weight,  once  a  charge  weight  of  200  pounds  is  exceeded. 
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This  is  a  data  report  which  presents  the  results  of  the  Navel  Ordnance  Laboratory 
(NGL)  blast  measurements  on  the  recently  completed  AN/FO  tests*  AN/FO  is  an 
explosive  mixture  of  ammonium  nitrate  and  fuel  oil.  It  is  being  developed  os  a 
TNT  replacement  for  large  scale  nuclear  sirblast  simulation. 

The  tests  were  conducted  in  cooperation  with  the  Defence  Research  Establishment 
Suffield,  at  Ralston,  Alberta,  Canada,  during  August  1969.  Several  other  U.  S. 
agencies  also  participated  in  this  program. 

This  effort  was  funded  by  the  Defense  Atomic  Support  Agency  through  the  Naval 
Strips  Engineering  Center  of  the  Naval  Ship  Systems  Command,  The  work  was  performed 
under  DASA  Subtask  NA  007-04,  Taak  NOL- 194 . 

The  use  of  company  names  throughout  this  report  is  for  technical  information 
purposes  only.  No  endorsement  or  criticism  is  intended. 


GEORGE  G.  BALL 
Captain,  USN 
Commander 

J.  PSi’KS 
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NOMENCLATURE 

A,  B,  C,  D,  E  =  coefficients  of  gage  calibration  fit  (Equation  (B-l)) 

Af  =  frequency  deviation  ,  (Hz) 

I  =  positive  impulse,  (psi-asec) 

PQ  *  atmospheric  pressure,  (psi) 

P  =  peak  side-on  overpressure,  (psi) 

p  a  instantaneous  overpressure,  (psi) 

R  =  distance,  (feet) 

Tq  a  atmospheric  temperature,  (°R) 

t  a  instantaneous  time,  (msec) 

TOA  a  shock  time  of  arrival,  (msec) 

W  a  charge  weight,  (pounds) 

0  a  time  decay  coefficient,  (Equation  (B-2)) 

a  a  pressure  decay  coefficient  (Equation  (£-3)) 

A  •  scaled  distance,  (feet /pound‘d) 

t  «  positive  duration,  (msec) 

SUBSCRIPTS  AND  SUPERSCRIPTS 

0  »  atmospheric  conditions 

1  «  standard  conditions,  (1^.?  psi,  519°R) 

2  »  test  site  conditions 

1  a  parameter  cube  root  scaled 
m  «  extrapolated  parameter 
ss  -  secondary  shock 
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X.  introduction 

1.1  Background 

There  is  a  continuing  requirement  for  the  development  of  a  large  ac&le 
field  source  for  generating  an  airblast  environment  that  simulates  that  of  a  nuclear 
explosion.  The  need  arises  from  programs  designed  to  test  the  vulnerability  and 
blast  hardness  of  military  hardware  and  strategic  structures.  For  example,  aa  part 
of  its  ship  blast-hardening  program  the  Navy  has  exposed  special  structures  e.g,, 
radar  masts,  deck  houses,  and  fully  operational  ships  to  particular  blast 
environments . 

Large  scale  simulation  techniques  with  charge  weights  up  to  500  tons,  usually 
involve  the  use  of  chemical  explosives.  The  explosive  for  such  a  source  should 
be  inexpensive,  easily  handled,  and  safe to  use  In  the  past,  raultiton  simulant 
charges  have  been  constructed  from  cast  TNT  blocks.  Some  work  hes  been  done  using 
balloons  filled  with  detonable  gases  (ref.  (l)). 

IXiring  19^7,  NOL  proposed  the  use  of  ammonium  nitrate/fuel  oil  (Ail/FO),  a 
commercial  blasting  agent,  as  a  replacement  for  the  TNT  used  in  large  blast  trials. 
Among  the  expected  advantages  Of  the  AN/PQ  system  were: 

1.  Increased  economy.  AJi/FQ  costa  about  seven  cents  per  pound  in  place  et 
ground  zero  compared  to  about  11$  to  $1.00  per  pound  for  TNT  in  place- -depending 
upon  whether  reclaimed  or  new  TNT  is  used. 

2.  The  availability  of  Afi/FO  in  large  quantities  and  at  convenient  locations; 
the  TNT  supply  is  limited  and  is  greatly  affected  by  munitions  requirements. 

3-  Fewer  Host  anomalies  (c.  g. ,  jetting  and  assymetrical  blast  fronts)  would 
be  expected  from  a  homogeneous  charge  in  contrast  to  the  block-built  TNT  charges , 

k .  Improved  safety.  The  nwaoniuia  nitrate  and  fuel  oil  are  non -explosive 
until  mixed.  Fully  mechanized  mixing  and  delivery  systems  are  already  developed 
for  charge  preparation  at  ground  zero. 

After  a  pioneering  effort  conducted  by  KOI  in  Nevada  during  i960  which  establishe 
the  detonability  of  unconfined  Aii/FO  and  which  determined  its  airblast  character¬ 
istics  (ref.  (2)),  the  way  was  paved  for  0  larger  scale  study.  In  the  fall  of  1^68 
N'OL  proposed  e  progr  sa  which  culminated  in  these  Aii/FO  trials  of  August  1969* 

1.2  Objectives 

The  general  objective  of  this  program  was  to  demonstrate  the  feasibility 
of  vising  Ah'/FO  as  the  explosion  source  in  the  Department  of  Defense's  nueleer 
airblast  vulnerability  and  hardening  program  This  objective  was  an  outgrowth  of 


1 


hoi®  70-32 


the  original  intent  to  satisfy  the  more  limited  blast  requirements  for  the  Navy's 
airblast-ship  hardening  program. 

A  number  of  specific  objectives  were  investigated  by  NOL  and  other  participating 

agencies  during  the  course  of  this  effort.  The  major  objective  was  to  determine 

\ 

the  alrblast  characteristics  of  AN/ TO.  Since  some  teat  structures  were  available 
from  earlier  tests  at  the  site,  a  secondary  objective  was  to  measure  their  blast 
response . 

The  primary  objectives  were: 

1.  To  verify  the  detonability,  scaling,  and  reproducibility  of  AN/FQ  for  charge 
weights  up  to  100  tons. 

2.  To  extend  existing  AN/FO  alrblast  pressure- time - distance  data  by  including 
measurements  from  close  to  the  charge  surface  out  to  the  1  pai  level. 

3*  To  study  fireball  growth  and  observe  blast  anomalies. 

4.  To  study  the  engineering  aspects  of  preparing  and  firing  bagged  and  bulk 
AN/TO  charges. 

5.  To  compare  the  alrblast  performance  of  AN/PO  with  TNT. 

6.  To  compare  the  blast  characteristics  of  bagged  and  bulk  20  ton  AN/FO 
hemispheres . 

7.  To  determine  the  temperature  stability  of  the  AN/ TO  explosive  prior  to  firing 
for  charge  weights  up  to  100  tons. 

3.  To  study  the  cratering  of  20  and  100  ton  AN/TO  hemispheres. 

n:ti  secondary  objectives  were: 

...  To  bias'  load  a  full  scale  frame  house  at  the  1.5  psi  level  from  the 
100  ton  AN/FO  test  (Event  III). 

2.  To  blast  load  an  underground  model  silo  and  buried  rock  inclusions  on 
Event  III. 

This  report  deals  with  NOL*  §  efforts  on  primary  objectives  1,  2,  l),  6,  and  7* 

The  other  agencies  participating  on  these  trials  will  report  their  project  results 
separately. 

1.3  Exp  crime  i]  tel  Program 

The  experimental  phase  of  this  program  was  carried  out  with  the  cooperation 
of  the  Defence  Research  Establishment  Suffield  (DRES),  at  Ralston,  Alberta,  Canada 
(ref.  (3))<  In  addition  to  the  field  support  they  provided,  ISES  made  shock  ticse- 
of -arrival  and  crater  measurements  as  well  as  high  speed  photographic  observations 
on  all  three  AN/FO  tests.  These  data  already  have  been  reported  in  reference  (4) 
and  will  be  Included  in  a  comprehensive  DABk  report  to  be  prepared  by  NOL  covering 
all  aspects  of  these  AH/TO  trials. 
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Other  U.  S.  Agencies  participating  in  these  trials  included  the  Ballistic 
Research  Laboratories  (BRL),  Naval  Civil  Engineering  Laboratory  (NCEL),  Naval 
Weapons  Center  (NWC)  and  the  U.  S.  Geological  Survey  (flSOS).  BRL  provided  pressure 
measurements  in  the  predicted  3000  to  30  psi  overpressure  range.  NCEL  made  studies 
of  the  response  of  rock  inclusions  and  made  body  motion  observations  of  a  model 
silo  exposed  to  the  blast  loading  of  the  100  ton  AN/PO  te3t.  NWC  made  observations 
on  a  two  story  frame  house  exposed  at  the  1,5  psi  level  from  the  100  ton  AN/FO  test. 
The  Geological  Survey  and  IKES  made  crater  studies,  NOL  was  responsible  for  the 
explosives  phase  of  these  trials  and  made  pressure  measurements  from  200  down  to  1  pi 
The  three  tests  conducted  during  the  Ali/iO  trials  at  DRES  during  August  1969 
were  as  follows : 

EVEN?  I  -  20  ton  ATI /F0  hemisplnere,  bagged.  Detonated  lb  Aug  1369. 

EVENT  II  -  20  ton  AIi/FO  hemisphere,  bulk  in  fiberglaa  shell. 

Detonated  21  Aug  19^9 • 

EVENT  III  -  100  ton  AR/FO  hemisphere,  bulk  in  fiborgles  shell. 

Detonated  23  Aug  1969, 

2.  Experiment  and  Procedtu*es 

2.1  Test  Site  and  field  Arrangement 

The  trials  described  in  this  report  were  conducted  at  the  Watching  Hill 
Blast  Range  of  the  Defence  Research  Establishment,  Sui'fleld,  at  Ralston.  Alberta. 
Canada.  Figure  1  shows  the  range  area  of  IKES. 

This  site  was  selected  for  reasons  of  logistics  and  because  it  enabled  direct 
comparisons  to  be  made  with  the  earlier  work  on  detonations  of  multiton  TF5* 
hemispheres  (ref.  (5)  sjtd  (6)).  The  physical  characteristics  of  this  site  have  been 
fully  described  in  the  Operation  Prairie  Flat  Operations  plan  (ref.  (?))• 

The  general  layout  of  the  ground  seroes,  liOL  cable  lines,  bunkers  and  camera 
positions  for  all  3  events  is  illustrated  in  Figure  2.  In  order  to  accommodate 
the  secondary  objectives  of  these  trials,  the  ground  zero  for  Event  III  ( 100  tons 
of  AN/yO)  was  selected  so  that  the  model  siloyinstruicentad  on  an  earlier  REST  test, 
could  be  blast  loaded  again.  The  frame  house  1700  feet  northeast  of  GZ  III  was 
therefore  exposed  at  the  1.5  psi  level.  This  house  was  repaired  after  being  exposed 
at  the  1  psi  level  on  Operation  Prairie  Flat.  The  ground  zero  for  Events  I  and  II 
were  placed  along  a  NE-SW  line  with  OZII  300  feet  from  GZIII,  and  HZ’®  I  and  II 
being  l60  feet  apart. 

KOL  had  0  gage  stations  on  each  event,  with  2  pressure  gages  at  each  station. 
(This  instrumentation  is  described  in  Appendix  A  of  this  report),  for  Events  I  and  I 
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tiie  NOL  gage  line  was  perpendicular  to  the  line  between  ground  zeroes.  This 
arrangement  permitted  an  easy  reorientation  of  the  gages  at  each  station  without 
establishing  new  stations.  The  200,  100,  50  and  20  psi  gage  stations  were  baffled 
flush  with  the  ground.  The  10,  5,  2  and  1  psi  gage  stations  were  above  the  surface 
and  used  disc  type  baffles.  Photos  of  both  types  of  gag*  stations  are  presented 
in  Figure  3* 

All  gage  cables  were  placed  in  a  trench  12-lQ  inches  deep.  The  cable  trenches 
connected  each  gage  station  to  a  common  NOL  cable  trench  which  ran  to  the  NOL 
instrumentation  trailer  some  3000  feet  distant. 

Most  of  the  ammonium  nitrate  for  the  AN/PO  was  delivered  to  the  Suffield, 
Alberta,  Canadian  Pacific  Reilroad  (CFR)  aiding  in  a  10  ton  hopper  car.  The 
siding  was  about  35  miles  from  the  test  site.  The  remainder  of  the  AN  was 
trucked  directly  to  the  range  in  22  ton  capacity  TRIMAC  tanker  trucks  from  the 
supplier  in  Calgary,  Alberta-^  Further  details  on  charge  construction  are  provided 
in  Section  2.2  of  this  report. 


2 . 2  Explosives  and  Charge  Construction 
2.2.1  AN/K?  Main  Charge 

The  main  charge  for  these  trials  was  a  9 Vo  by  weight  AN/ H)  mixture  £.?4$ 
ammonium  nitrate  and  6$  fuel  oilj .  The  AN  itself  was  a  commercial  fertilizer 'and 
was  basically  the  same  type  of  grills2  used  in  our  1968  Nevada  tests  (ref.  (2)}. 
The  FO  was  summer  grade  No.  2  diesel  fuel.  A  rod  di  e  was  added  to  the  fuel  oil 
to  pemit  •  ready  visual  check  on  the  AN/ FO  mixing  proportions. 

The  20  tons  of  bagged  AN/PO  for  Event  1  were  prepared  at  the  01  area .  Tim 
AN  was  transported  from  the  ?0  ton  hopper  cor  at  the  Suffield  CPR  siding  by  the 
AN/PO  mixing  truck.  The  truck  had  a  capacity  of  about  7  tons  of  mixed  product 

•5 

(Ali/FO).  Thus,  three  loads  Wei's  required  per  20  ton  event.  A  bagging  unitJ  was 


explosives  contractor  to  NOL  was  Ace  Explosives  Ltd  of  Calgary,  Alberta.  The 
AN  used  was  manufactured  by  Cociiuco  Ltd  also  of  Calgary,  Alberta. 


^Frills  are  poruu®,  spherical  particles.  They  arc 
in  a  prilling  tower  and  are  much  like  lead  shot  in 
density  o?  about  1.4  g»/cc  compared  to  the  crystal 

"The  bagging  unit  was  designed  and  built  by  Hr.  0. 


formed  by  dropping  molten  AN 
size  and  shape.  They  have  a 
density  of  AN  which  is  1.725  gm/cc. 

R.  Rintoul  of  Ace  Explosives,  Ltd. 


4 


NOI/ER  70-32 


located  at  the  site  during  Event  I  charge  placement.  The  arrangement  of  the  mixer 
truck  and  bagging  unit  at  the  GZ  location  is  shown  in  Figure  4. 

A  total  of  800 -50  pound  bags  was  prepared  for  this  first  80  ton  charge.  Six  hun 
dred  and  fifty  50- lb  bags  were  used  in  the  layered  arrangement  illustrated  in  Figure 
Loose  AN/PO  from  150  of  the  bags  was  poured  into  the  spaces  between  bags  to  form 
a  charge  with  uniform  density,  i.  e . ,  no  airspaces.  The  bag  dimensions  were 
21  x  13-5  x  5*o  inches. 

The  Ak/F0  mixer  truck  used  a  system  of  augers  to  feed  the  AN  from  the  bins  to 
the  fuel  metering  point.  At  this  point  the  red-dyed  fuel  oil  was  mixed  with  the 
AN.  The  mixed  AN/FO  was  then  fed  through  a  vertical  auger  and  out  a  swinging 
horizontal  auger  to  place  the  AN /Ft  where  it  was  needed,  i.e.,  into  the  bagging 
unit  for  Event  I  and  into  the  charge  cases  for  Events  II  and  III. 

Changes  in  the  fuel  oil  content  of  the  AN/FO  were  detected  very  quickly  by 
visual  means  because  of  the  red-dyed  fuel  oil  used.  The  fuel  oil  content  was  also 
monitored  quantitatively  throughout  the  explosives  placement  operation  by  chemical 
analysis  (ref.  (9))-  Table  1  contains  data  on  AN/FO  charge  dimensions,  weight, 
density  mid  fuel  content  for  all  three  events. 

Events  II  and  III  were  charges  of  bulk  AN/FO  placed  into  thin  walled  fiberglas  / 
polyester  resin  containers.  The  fiberglas/ polyester  resin  shell1  was  made  of  section 
having  full  compound  spherical  Curvature . 

The  80  ton  size  container  fur  Event  II  was  14.0  feet  in 

,  p 

base  diameter  and  was  made  of  11  section**  each  3 Ah  inch  thick  .  The  100  ton  sise 
container  was  24.2  feet  in  base  diameter  and  had  22  sections, each  1/4  inch  thick.  . 

Tiie  sections  were  joined  together  by  nylc  bolts  and  epoxy  resin  adhesive. 

To  fill  the  fiberglas  shells  the  sixer  truck  was  backed  up  to  the  container 
for  AN/FO  placement.  Hie  time  required  for  each  mixing  and  placement  cycle  (i.e., 
each  7  tons  of  AN/EO)  was  about  three  hours .  To  reduce  the  loading  time 
on  the  lOo  ton  AS/KO  charge  cost  of  the  AN  was  trucked  from  the  Coainco  fertiliser 
plant-  in  Calgary  to  the  GZ  via  22  ton  capacity  TRIKAC  tanker  trucks.  The  AN  was 

1  The  sections  for  the  two  containers  were  manufactured  by  i\ ogay  Itodel3  of 
Betheada,  Mai 7  land . 

2 

Before  the  field  operation,  vc  conducted  high  speed  camera  testa  on  samples 
of  the  shell  naterlal  to  determine  their  behavior  when  in  contact  with  detonating 
explosives.  The  high  speed  photographs  indicated  break-up  of  the  fiberglas 
within  a  few  inches  of  the  charge. 
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fed  into  the  mixer  truck  and  the  AN/FO  into  the  container  in  a  continuous  operation. 
Using  this  loading  system,  about  23  tons  of  AN/FO  were  mixed  and  placed  in  a 
four  hour  period.  A  photograph  illustrating  the  arrangement  of  the  TRIMAC  tanker 
and  AN/FO  mixer  trucks  at  the  100  ton  G2  is  presented  in  Figure  6.  The  completed 
20  ton  audio 0  ton  bulk  AN/FO  charges  are  shown  in  Figures  7  and  8. 

2.2.2  Booster  and  Primacord  Initiation  Method 

The  hemispherical  boosters  used  for  all  three  events  were  prepared  by 
the  U.  S.  Naval  Ammunition  Depot,  Hawthorne,  Nevada  (ref.  (8)).  The  boosters 
were  a  nominal  250  pounds  each  total  weight  and  consisted  of  a  16  pound  hemispherical 
50/50  pentolite  primer  with  about  23^  pounds  of  TNT  cast  over  it. 

NOL  developed  a  primacord  initiation  method  (ref.  (2))  which  was  used  for 
each  event.  In  this  method  a  strand  of  100  grains  per  foot  primacord  is  placed 
in  a  shallow,  radial  trench  beneath  the  charge,  leading  from  the  GZ  to  beyond  the 
outer  edge  of  the  AN/FO  charge.  The  C-Z  end  of  the  primacord  is  fed  through  a 
radial  hole  in  the  booster  and  a  small  knot  is  tied  at  the  top  to  secure  it. 

This  method  greatly  simplified  the  arming  procedure,  as  the  electric  detonator  is 
simply  attached  to  the  other  end  of  the  primacord  still  exposed  after  the  charge 
has  been  completed.  The  explosive  train  is:  electric  detonator  —♦.primacord  — *. 
pentolite  primer  -*.TNT  booster  —*•  main  charge  (AN/FO).  Tula  la  all  illustrated 
schematically  in  Figure  9 • 

2 . 3  NOL  Instrumentation 

Airblast  pres suit:  histories  were  measured  with  variable  reluctance  transducers 
and  recorded  on  magnetic  tape  recorders.  The  temperature  within  the  AN/FO 
charge  was  monitored  with  thermistors  on  each  event.  The  instrumentation  system 
is  described  la  detail  iu  Appendix  A. 

Before  each  event,  both  the  pressure  gages  and  the  thermistors  were  statically 
calibrated.  The  calibration  of  the  thermistor  took  into  account  the  resistance 
of  ‘the  cable  between  ground  sero  and  the  Instrumentation  trailer. 

2.k  Data  Analysis  Procedures 

The  pressure-time  records  were  digitized  and  then  analysed  using  techniques 
which  are  described  in  detail  in  Appendix  B.  Hie  parameter*  computed  include 
peak  pressure,  positive  duration  and  positive  impulse.  Extrapolations  to  peak 
pressure  end  positive  duration  were  Bade  using  techniques  described  by  Ethridge 
(ref. (10).  These  extrapolation  techniques  were  used  to  take  into  account  both 
the  finite  rise-tiae  c-f  the  observed  signal  incurred  because  of  instrumentation 
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system  limitations  and  also  any  early -time  gage  malfunctions.  Gage  malfunctions 
were  noted  on  several  signals  from  each  event- -namely,  at  those  stations  in  the 
50  psi  region  and  above.  These  malfunctions  were  manifested  by  a  loss  of  PM 
carrier  amplitude  for  several  milliseconds  upon  the  arrival  of  the  airblast  wave. 
This  loss  of  carrier  exhibited  itself  as  spurious  peaks  on  the  discriminated  signal. 
To  handle  this,  only  that  portion  of  the  record  that  occurred  after  the  gage 
resumed  normal  operation  was  used  in  the  computations. 

A  comparison  between  the  procedure  described  in  Appendix  B  for  determining 
impulse  and  a  direct  measurement  with  a  planimeter  on  several  pressure-time  records 
showed  excellent  agreement  (within  a  few  percent). 

Time  of  arrival  data  were  measured  directly  from  the  tape  recordings,  using 
an  electronic  counter  operating  in  the  time  interval  mode. 

3.  Airblast  Results 

3.1  The  Data  and  Scaling  Procedures 

The  unsealed  data  obtained  on  all  three  events  are  presented  in  Tables  2,  3, 
and  k.  The  data  presented  in  these  tables  are  for  the  ambient  conditions  at 
DRES  as  shown  in  Table  5.  To  make  for  a  useful  comparison  with  previously 
published  TNT  data  (ref.  (5)  and  (6)),  the  AN/PO  data  were  cute  root  and  Sachs' 
Scaled  (ref.  (ll))  to  standard  sea- level  conditions  of  pressure  and  temperature. 

To  do  this  scaling,  the  following  equations  were  used: 


For  Pressure: 


(1) 


j 

For  Distance: 


(2) 


For  Times: 


(3) 
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And  for  Impulse: 


The  scaling  factors  for  all  three  events  are  presented  in  Table  5. 

The  data  fbr  each  station  of  Tables  2,  3  and  4  were  averaged  and  the  scaling 
equations  (Equations  1  through  4)  vert  applied.  The  resulting  scaled  averaged 
data  are  presented  in  Tables  6,  7  and  S  for  Events  X>  II  and  III  respectively. 

All  tabulated  data  are  given  to  three  significant  figures. 

The  peak  pressure  versus  scaled  distance  data  (p  vs  A  )  for  all  three  events 
are  shown  graphically  in  Figure  10.  The  TNT  standard  curve  (ref.  (5))  is  also 
plotted  in  this  figure  to  enable  the  making  of  direct  comparisons  between  AN/FQ 
and  TNT. 

A  5th  degree  polynomial  was  fitted  to  a  composite  of  all  of  the  pressure  (P1t1)- 
scaled  distance  ( \)  data  of  Tables  6,  7  and  8.  Tire  form  of  the  equation  was; 

In  .  10.4781  -  9.01448  (In  \)  +  5-55124  (In  \)2  -  2.33879  (in  A)3  + 
.514723  (In  A)4  -  .0447655  (In  A)?  (5) 

Thin  equation  is  ’/olid  over  the  1  to  200  psi  region  and  is  represented  by 
the  solid  line  in  Figure  10. 

Figure  11  is  a  plot  of  tiro  scaled  positive  duration  and  scaled  distance 
( ^  vs  A)  data.  The  scaled  positive  impulse  --  scaled  distance  data  (i^  vs  A) 
ore  plotted  in  Figure  12.  a  TNT  standard  curve  from  ref.  (6)  is  also  plotted  in  Fig.  12. 

3 • 2  Secondary  Shock  Measurements 

A  late  secondary  shock  wave  was  measured  on  the  20  psi  and  below  pressure 
records  on  all  three  events.  This  distinct  secondary  shock  occurred  during  the 
negative  phase  portion  of  the  pressure-time  curves  (see  Appendix  C). 

Information  on  the  secondary  shock  is  seldom  reported,  although  a  review  of 
earlier  work  reveal a  its  presence  on  a  majority  of  the  data  records  ( e .  g . ,  ref.  (12)). 
Because  of  its  hydrodynamic  interest  and  its  possible  significance  for  response 
test  applications,  secondary  shock  information  is  included  in  this  report,  lire 
wa-i-e  shapes  can  be  observed  on  the  records  reproduced  in  Appendix  C.  The 
secondary  shock  amplitudes  and  times  of  arrival  with  respect  to  the  detonation  zero 
pulse  are  presented  in  Tables  2,  3»  6,  7  and  8.  figure  13  is  a  plot  of  all  the 
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time  of  arrival  data  (main  shock  and  secondary  shock)  as  a  function  of  scaled  distant 
The  scaled  peak  secondary  shock  pressure  is  plotted  versus  scaled  distance  for  all 
3  events  in  Figure  14. 

3.3  Temperature  Measurements  in  the  AN/FO 

The  temperature  within  the  AN/FO  explosive  was  monitored  with  a  thermistor. 

For  each  event,  the  thermistor  was  placed  at  the  center  of  mass  of  the  hemisphere. 

It  was  felt  that  this  location  would  be  warmest  if  any  self  heating  was  to  take 
place  in  the  AN/FO.  For  the  20  ton  charges  the  thermistor  was  located  3  feet  above 
the  base  of  the  hemisphere.  Similarly,  for  the  100  ton  charge  the  thermistor  was 
located  at  5  feet  above  the  base  of  the  hemisphere.  In  addition,  as  a  second 
check,  another  thermistor  was  placed  at  the  top  of  the  booster  on  Event  III. 

The  temperature  change  in  the  AN/FO  mass  was  of  negligible  magnitude.  The  tem¬ 
perature-time  data  are  plotted  in  Figures  1 5,  l6  and  17  for  Events  I,  II  and  III 
respectively.  Note  the  cooling  down  of  the  explosive  at  times.  At  no  time  did 
the  recorded  temperatures  exceed  o7“F.  The  initial  temperature  of  the  AN  in  the 
hopper  car  and  in  the  tanker  trucks  was  about  88°F.  The  ambient  air  temperature 
reached  a  maximum  of  105°F  during  the  loading  of  the  100  ton  AN/FO  charge  for 
Event  HI. 

4 .  Discussion  and  Conclusions 
4 . 1  General 

All  three  charges  detonated  properly  and  high  order.  This  is  evidenced  by 
the  following  observations: 

a.  The  pressure-time  records  show  the  familiar*  and  classical  wave  shapes 
(see  Appendix  C). 

b.  The  results  of  the  20  ton  shots  scale  well  with  the  100  ton  data 
(see  Figures  10-12).  The  extent  of  data  scatter  is  of  the  same  order  as  for  TOT 
fired  under  the  same  conditions  and  is  attributed  to  a  large  extent  to  the 
accuracy  of  the  instrumentation  and  to  vagaries  of  field  operation. 

c.  The  results  of  these  large  scale  urbu.u  scale  well  with  the  earlier 
260  lb  to  4000  lb  AH/FO  results  (ref.  (2)). 

From  these  observations  it  can  be  deduced  that: 

a.  Over  the  pressure  range  measured  by  this  project,  there  is  no  signifi¬ 
cant  difference  between  the  tagged  and  bulk- loaded  AH/FO  charge®. 
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b.  The  fuel-oil  does  not  settle  out  of  the  AN/FO  mixture;  if  this 
lad  occurred,  it  could  he  expected  that  the  tvo  20  ton  shots  would  have  given 
different  results  and  would  not  have  scaled  to  the  100  ton  and  small  size  charge 
results,  (indeed,  visual  observation  during  charge  preparations  and  prior  to 
firing  time  did  not  show  any  evidence  of  oil  seepage.) 

c .  From  Figures  10-12, it  is  evident  that  there  is  no  significant  difference 
between  the  pressure-distance  characteristics  of  AN/FO  and  TUT. 

A  single  frame  from  one  of  the  Canadian  high  speed  camera  films  on  Event  X 
is  shown  in  Figure  18 .  This  photograph  shows  the  smoothness  and  symmetry  of  the 
shock,  wave  (at  42.9  milliseconds  after  detonation)  produced  by  the  20  ton  bagged 
AN/FO  hemisphere . 

Forzei  (ref.  (15))  and  Lehto  (ref.  (16))  have  independently  made  calculations 
of  the  pressure-distance  characteristics  of  the  AN/FO  system.  Both  sets  of 
analyses  show  good  agreement  with  our  experimental  results  in  the  1-200  psi  region. 

The  second.'try  ehock  (described  in  section  3-2)  which  oceurs  near  the  minimum 
of  the  negative  phase  in  Large  explosion  trials  deserves  further  attention.  For 
structural  elements  exposed  at  low  pressures,  the  secondary  shock  could  be  very 
importafct.  This  is  because,  as  the  m?iu  shock  propagates  and  decays,  the  ratio 
of  its  amplitude  to  the  secondary  shock  amplitude  tends  to  approach  a  value  of  one. 

4.2  Equivalent  Weight  of  AN/FO 

Long  standard  NOL  procedures  for  evaluating  th_  peak  pressure  and  impulse  TNT 
equivalent  weights  (EW^  and  KWj.)  (ref.  (13)  and  (l4);  were  used  on  the  preeert  data. 

A  composite  of  the  data  presented  in  Tablet  6,  7  and  8  was  used  in  this  analyola . 

The  pressure  and  impulse  versus  distance  ourvos  for  any  set  of  tee*  and  standard 
explosives  are  seldom  parallel.  Thus,  the  slngl  Tkiue  of  average  equivalent 
weight  usually  given  for  a  test  explosive  may  be  misleading  because  it  cannot 
indicate  hov  it  varies  as  a  function  of  side-on  overpressure.  To  illustrate  this 
functional  variation,  the  EW  add  EWT  for  Alt/ 10  are  plotted  as  a  function  of 
pressure  in  Figure  1$. 

It  is  interesting  to  note  that  although  EVL.  varies  from  0.77  to  1.17, this 
magnitude  of  variation  in  yield  is  hardly  suspected  when  viewing  th..  pressure - 
distance  curves  of  Figure  10*  Figure  10  shows  a  scatter  of  data  around  either  the 
TNT  or  AH/JO  curves  nc  greater  than  that  usually  observed  on  large  scale  field 
trials. 
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Equivalent  weight  determinations  are  an  exceedingly  sensitive  measure  of  the 
merits  of  tone  explosive  compared  to  another;  for  some  applications  it  may  he  too 
sensitive  and  hence  of  little  practical  significance.  And  of  even  lesser  practical 
significance  may  he  the  averaged,  single-value  equivalent  weight  —  particularly  if 
the  average  is  taken  over  a  large  pressure  range.  The  concept  of  averaged,  single- 
valued  equivalent  weights  is  so  rooted  in  the  explosives  field,  however,  that 
although  it  is  with  trepidation*  we  present  these  values.  The  user  must  observe 
the  pressure  range  over  which  the  averages  are  taken  and  he  aware  of  the  limitations 
of  these  average  values. 

*t 

The  average  EWp  for  AN/FO  over  the  1  to  200  psi  range  is  0.94  *  .06  relative 
to  TNT  using  a  logarithmic  weighting  method.  Similarly,  the  average  EW^.  for  AN/FO 
is  0.71  t  .05.  Using  logarithmic  averaging  over  the  1-30  psi  range  (the  data  range  < 
earlier  AN/FO  work  (ref.  (2)),  the  average  EWp  is  0.86  ±  .03.  In  reference  (2)  we 
used  a  linear  weighting  method  when  we  averaged  the  equivalent  weights  and  arrived 
at  a  figure  of  0.82  for  the  EWp.  The  linear  weighting  gives  greater  emphasis  to 
the  EMp  at  the  higher  pressure  levels.  Using  the  present  logarithmic  method  on  the 
phase  I  AN/FO  Data  of  reference  (2),  an  average  SWp  of  O.87  is  obtained  over  the 
1-30  psi  range. 

4.3  Thermal  Stability 

Hie  AH/PO  temperature  data, as  presented  in  Figures  15,  16  and  17,  indicate 

that  massive  AN/FO  has  good  thermal  stability.  In  the  Event  I  data  there  Is  a 
definite  cooling  trend.  In  the  case  of  the  Event  II  and  III  data  there  is  some 
evidence  of  a  very  slight  general  warming  trend  among  the  cooling  and  warming 
cycles.  It  is  on  the  order  of  about  IF0  per  day  (the  measurements  are  accurate  to 
within  tlF°).  At  this  point  It  would  appear  to  be  not  self  heating  of  the  AN/FO 
but  rather  external  heating  from  the  sun .  During  the  loading  of  the  Event  III 
charge  the  outside  temperature  reached  about  105°F.  AN/FO  Is  a  good  insulator. 

The  initial  temperature  of  the  AW  was  about  Q8°F. 

It  can  be  concluded  that  the  AN/FO  did  not  exhibit  any  self  heating.  No  self 
heating  is  expected  for  50°  bon  AN/FO  chargee. 

4.4  Concluding  Statement. 

The  results  of  these  AN/FO  trials  in  conjunction  with  NOL'a  earlier  work 
on  AN/PO  have  built  up  a  now  formidable  portfolio  of  data  on  the  airblast  performance 

1  - * - 

The  standard  deviation  of  the  mean. 

2 

Averaging  the  values  of  equivalent  weight  at  logarithmic  pressure  intervals 
(i.e.,  1,  2,  5,  10,  20,  50,  100  and  200  psi). 
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of  this  explosive .  AN/FO  offers  advantages  of  economy,  safety,  ease  of  handling, 
availability  and  reproducibility  over  TNT>  slurried  explosives,  or  any  other  system 
presently  in  use  for  large  scale  simulation  of  nuclear  air  blast. 

With  these  advantages,  AN/FO  can  be  seriously  considered  as  a  candidate  explosive 
for  future  large  chemical  explosive  trials. 
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FIG.  1  THE  DEFENCE  RESEARCH  ESTABLISHMENT,  SUFFIELD  RALSTON 
ALBERTA,  CANADA 
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(a)  FLUSH  BAFFLE.  NOL  GAGE  ON  LEFT. 
BRL  GAGE  ON  RIGHT. 


(b)  STANDOFF  BAFFLE 

FIG.  3  AIRBLAST  GAGE  MOUNTS 


EVENT  L  MIXING  T.iUCK  IS  ON  RIGHT;  SAGGING 


COMPLETED  CHARGE  FOR  EVENT  I 
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TANKER  TRUCK  ON  l£:T  FEEDING  AN  INTO  AN/FO  MIXER  TRUCK.  AN/FO  ENTERING  CONTAINER  FOR  EVENT®. 
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FIG.  9  SCHEMATIC  ARRANGEMENT  OF  THE  AN/FO  CHARGES  OF  EVENTS  H  AND  HI 


OVERPRESSURE 


FIG.  10  PEAK  PRESSURE  VERSUS  SCALED  DISTANCE  FOR  AfVFO  AT  ORES.  SCALED 
TO  SEA  LEVEL  CONDITIONS. 


FIG.  11  SCALED  POSITIVE  DURATION  VERSUS  SCALED  DISTANCE  FOR  AN/i'O  AT  ORES. 
SCALED  TO  SEA  LEVEL  CONDITIONS. 
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V  A 

V 
□ 

°\\ 

O  EVENT  I  20  TON  BAGGED 

A  EVENT  H  20  TON  BULK 

Q  EVENT  nr  100  TON  BULK 

-  AN/FO  COMPOSITE 

- TNT  STANDARD  (REF  6) 

\ 

_ 1 _ 1 _ 1 _ 1 _ L _ LJ— L 

1  JO  1( 

1/3 

SCALED  DISTANCE,  k  (FT/L8  '  ) 


FIG.  12  SCALED  IMPULSE  VERSUS  SCALED  DISTANCE  FOR  AN/FO  AT  DRES.  SCALED 
TO  SEA  LEVEL  CONDITIONS. 


SCALED  TIME  OF  ARRIVAL  (MSEC/LB 
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FIG.  13  SCALED  TIMES  OF  ARRIVAL  VERSUS  SCALED  DISTANCE  FOR  AN/FO  AT  DRES 
SCALED  TO  SEA  LEVEL  CONDITIONS. 
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1  10 

SCALED  DISTANCE,  A(FT/LB1,/3) 


FIG.  14  PEAK  SECONDARY  SHOCK  PRESSURE  VERSUS  SCALED  DISTANCE.  AN/FO  TRIALS, 
ORES,  AUGUST  1969 


FIG.  15  TEMPERATURE-TIME  HISTORY  IN  Ahj/FO,  EVENT  I 


TEMPERATURE-TIME  HISTORY  IN  AN/FO,  EVENT  H 
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FIG.  19  PEAK  PRESSURE  AND  IMPULSE  EQUIVALENT  WEIGHT  VERSUS  PRESSURE 
AN/FO  AT  ORES,  AUGUST  1969 
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TABLE  3  TABLE  4 

NOL  AIR3LAST  MEASUREMENTS,  E'/ENT  2,  UNSCALED  DATA  NOL  A1RBLAST  MEASUREMENTS,  EVENT  HI,  UNSCALED  DATA 

W  =  37,600  LBS  W=  200,900  LBS 
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APPENDIX  A 
HOL  INSTRUMENTATION 

The  pressure  gages  used  In  these  tests  vere  variable  reluctance  transducers 
manufactured  by  Consolidated  Controls  Corporation.  These  are  frequency  nodulated 
(1W)  gages  which  operate  in  the  standard  IRIG  13  and  IRIG  14  frequency  bands, 

14.5  has  and  22.0  kHz  respectively. 

The  gage  signals  were  transmitted  to  the  instruiaentation  trailer,  some  3000  fee 
from  the  u.  Z.  (s,  over  WDl/TT  field  telephone  wire.  The  signal  oablea  were 
terminated  by  United  Transformer  company  model  UTC  A-12  transformers  and  the 
signals  then  recorded  on  magnetic  tape  recorders.  Three  14 -track  recorders  were 
used:  1)  Amp ex  FP.  l800L,  2)  Consolidated  Electrodynamics  Corporation  VR  3300,  and 
3)  Sangamo  4700.  The  FM  signals  were  all  recorded  in  the  direct  record  ( amplitude 
modulated)  mode. 

Pressures  were  measured  at  eight  distances  on  each  event,  with  two  pressure 
transducers  at  each  distance.  A  time  zero  pulse,  provided  by  the  DRES  control 
bunker,  was  also  recorded  on  each  shot.  Thus,  17  channels  of  information  were 
recorded  on  each  event.  The  incoming  signals  were  divided  in  such  a  way  that 
any  two  of  the  three  recorders  contained  a  complete  set  of  records 

On  playback,  the  signals  were  played  through  a  tunable  discriminator  manufaetur 
by  Slectro-Meehanical  Research  Corporation  and  recorded  on  a  Midwestern  oscillograp 
The  oscillograph  records  of  pressure  versus  time  were  digitized  using  the  KOI 
Telereader  system. 

The  system  frequency  response  was  flat  from  D.C.  to  1  kHz.  the  gages  being 
the  response -limiting  element.  This  relatively  low  upper  frequency  response  was 
sufficient  for  the  long  duration  signals  expected  and  observed  on  these  trials. 

This  low  upper  frequency  response  manifests  itself  as  a  finite  rise-time  and  a 
reduction  in  the  apparent  peak  amplitude  of  the  observed  gage  signalE.  Using  zh  * 
extrapolation  techniques  described  in  Appendix  B,  the  observed  signals  are  extra¬ 
polated  back  to  zero  time  (shock  arrival).  This  procedure  corrects  for  the  upper 
frequency  limitations  of  the  system. 
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The  temperature  of  the  explosive  in  each  charge  was  monitored  by  thermistor*. 
A  General  Radio  Type  1650-A  Impedance  Bridge  va*  used  to  read  the  thermistor 
resistance.  The  accuracy  of  this  measuring  system  was  ±1F°. 
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APPENDIX  B 

DATA  ANALYSIS  PROCEDURES 


A  least  squares  fit  of  the  fora 

P  -  A(tf)  +  B(Af)2  +  C(tf)3  +  B(AfA  +  ,  (B-i) 


was  made  to  the  calibration  data  for  each  gage  for  each  event.  A  3rd,  hth  or  3th 
degree  polynomial  vas  chosen  for  each  set  of  calibration  data.  The  smallest  degree 
of  fit  for  optimum  accuracy  vas  selected  for  each  set  of  calibration  data. 

The  digitized  data  for  each  P-t  record,  along  with  the  coefficients  of  the 
gage  calibration  data  ( Equation (B-l )) was  analyzed  using  the  IBM  7090/709^  computer. 
The  calculational  methods  used  and  a  listing  of  the  computer  program  are  presented 
herewith  on  pages  B-3  to  B -6. 

Extrapolated  positive  duration  was  determined  by  fitting  an  equation  of  the 
form: 

t  .  (  (B-2) 

to  the  pressure-time  data  in  the  last  quarter  of  the  apparent  positive  phase.  The 
value  of  Tg^  is  the  extrapolated  positive  duration. 

Extrapolated  peak  pressure  was  determined  by  fitting  an  equation  of  the  form: 

p  Pgm®0*1  (B-3) 


to  the  pressure-time  data  in  the  first  half  of  the  apparent  positive  phase.  The 
value  of  Pg^  is  the  extrapolated  peak  pressure . 

Positive  Impulse  is  defined  by  the  equation; 

/*T 

I2  •  J  2  p(t)  dt.  (B-UA) 


In  these  calculations,  the  impulse  was  determined  in  two  parts .  Over  most  of  the 
positive  phase,  after  some  initial  time  interval  At,  the  impulse  vas  determined 
by  the  equation.  _ 


"at 


p(t)  dt, 


(b-4b) 


where  At  is  a  small  value  of  time,  which  accounts  for  both  the  rise-time  of  the 
ob a erred  signal  and  any  observed  early- time  gago  malfunctions.  Over  this  range 
(At  to  tg),  the  impulse  was  determined  by  the  use  of  the  trapezoid  rule  --  that  is, 

a  numerical  integration  of  the  pressure -time  data. 

B-J 
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The  Impulse  In  the  time  increment  (At)  between  shock  arrival  and  the  first 
* 

pressure  point  was  determined  in  the  following  way. 


P  -  , 

(B-3) 

/.At 

Al  ~j  p(t)  dt  t 

(B-4c) 

■o  /-At  . 

AI  aP2»  f  eat  dt  , 

Jo 

(B-4d) 

(2-416) 

This  impulse  increment  (Equation  (b-4i))  was  then  added  to  the  impulse  determined 
for  the  remainder  of  the  positive  phase  to  arrive  at  the  total  positive  Impulse 
(that  is  I_  a  I  +  AI). 
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COMMON  X(4,500)*A(10> »P ( 500 ) *LL ( 50 > »TL ( 500) .PT ( 500 ) tTPLOT ( 500) 
1PL! 500) .DUMMY (50)  *1X(2»  500)  »T  ( 500) »U( 500) *TI TLE(24)  *01. G! 500) 

JCASE  IS  THE  NUMBER  OF  RECORDS  BEING  PROCESSED 
READ ( 5 • 5000 ) JCASE 
DO  999  KIK*1 » JCASE 
READ( 5 *5100) (TITLE! I >.1*1.4) 

WRITE(6»5110)CTITLE(I).I*1.4) 

IDEG  IS  THE  DEGREE  OF  THE  POLYNOMIAL  USED  TO  FIT  THE  CALIBRATION 
DATA  FOR  THAT  GAGE  AND  SHOT 
READ(5.5000) IDEG 

XCAL  AND  YCAL  AR5  THE  SIZE  OF  THE  X  AND  Y  CALIBRATION  STEPS. 

READ! 5  *5140) XCAL. YCAL 

XCAL  IS  IN  MILLISECS  AND  YCAL  IS  IN  HERTZ. 

IXSCA  AND  IYSCA  ARE  TELEREADEX  CALIBRATIONS 
READ! 5  *5130 ) I XSCA  » I YSCA 
XSCA-IXSCA 
YSCA*I YSCA 

I TOA  IS  THE  SHOCK  TIME  OF  ARRIVAL  OBTAINED  FROM  THE  RECORD 
IDUM  IS  A  DUMMY  VARIABLE 
READ (5*5130) ITOA. IDUM 
TOA* I TOA 
DUMMY* IDUM 

THE  A! J )  ARE  THE  COEFFICIENTS  OF  THE  CALIBRATION  CURVE  FIT 
DO  10  Jsl.IDEG 
10  READ ( 5* 5150 ) A ( J ) 

M*0 

IX(l.L)  AND  IX ( 2 »L )  ARE  THE  POINTS  PUNCHED  BY  THE  TELEREADER 
SYSTEM 

DO  20  L-1.500 

READ ( 5.5130) IX! 1»L)*IX(2*L) 

X  ( 1 »L ) » I X ( 1 » L  ) 

X(2*L)-IX(2.L) 

I F ( X ( 1 *L ) .EO. 999999 . ) GO  TO  25 
M*M+ 1 

20  CONTINUE 
25  MM«M 

C 

WR I TE  <  6 .5200 ) MM 
XS1*A8S!XSCA) 

XS2  *ABS ( XCAL  > 

YSl-ABS(YSCA) 

YS2»ABS( YCAL) 

C 

DO  40  JJ»1.MM 

C  X ( 3  .  JJ )  IS  THE  TIME  CALCULATED  FOR  EACH  POINT 

X(3.JJ)*(X<1 »JJ)/XS1)*XS2 

C  X ( 4* JJ )  IS  THE  FREQUENCY  DEVIATION  CALCULATED  FOR  tACH  POINT 

X(4.JJ)*(X(2.JJ)/YS1)*YS2 
T( JJ ) »X ( 3. JJ ) 

U( JJ)*X(4.JJ) 

40  CONTINUE 
C 

WRI TEI6.5210) 

DO  90  K* 1 »MM 

C  P(K)  IS  THE  OVERPRESSURE  CALCULATED  FROM  EACH  FREQUENCY  DEVIATION 
P ( K ) *0. 

DO  80  KK»1 . IDEG 
P  !  K  )  «P  (  K )  4-A!  KK )  *X  (4  .K  )  * *KK 
80  CONTINUE 
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90  CONTINUE 
C 

SI1)>0« 

NNN-MM-1 

C  THIS  SECTION  CALCULATES  IMPULSE  BY  THE  TRAPEZOID  RULE 
00  200  LIU1.MM 
IN-LK-1 

DELT»ABS( T ( LK+1 )-T(  Lit ) ) 

IF  ( P(  LK. )  .EQ.O.  )G0  TO  150 
GO  TO  155 
150  G(  Lit)  *0. 

GO  TO  160 

155  G(LK)sG( IN)+«5*(P(LX)+P{ LK+1 ) ) *0ELT 
160  WRITE(6.522Q)T<LX).P<LX).G(LX> 

200  CONTINUE 

CALL  OURAT ( T.U.P.MM ) 

OURAT  DETERMINES  BOTH  THE  ACTUAL  CROSSING  TIME  OF  THE  SIGNAL  AND 
ALSO  THE  EXTRAPOLATED  DURATION 

CALL  PRES1 ( Dl.T  »P»MM) 

PRES1  CALCULATES  THE  EXTRAPOLATED  PEAtt  PRESSURE 

999  CONTINUE 
5000  FORMAT  a  15) 

5100  FORMAT ( 4A6 ) 

5110  FORMA T{ 1H1 *4A6 ) 

5130  FORMAT (117*1110) 

5140  FORMAT ( 2E1 0.4  > 

5150  FORMAT <  EX4 • 5 > 

5200  FORMAT ( 1H0  »  2HM« » 1 15 ) 

5210  FORMAT(1H0.44HTIME(MSEC)  PRESSUREl PSI )  IMPULSE ( PS I-MSEC >•✓/ > 

5220  FORMAT ( 3F10«4 ) 

STOP 

END 

J1BFTC  SDURA 

SUBROUTINE  DURAT  (T.U.P.MM) 

COMMON  X (4.500). A( 101 »P < 500 > »LL ( 50 ) » TL ( 500 ) *PT ( 500 ) .TPLOT(5Q0>  • 

1PL  ( 500 )  .DUMMY I  50)  *IX( 2*500)  »T <  500 > .U< 500 ) *T I TLt (24)  .Ui  .6(500) 

DO  30  NN-1.50 
LL(NN)«0 
30  CONTINUE 
X*  1 
Xl«0 

DO  430  1*1. MM 
IF  <  U( 1 ) >430.420.430 
420  LL(K)«I 
X»X+1 
X I »K I  +  1 
430  CONTINUE 

IF  <  LL (2 )-Xl 1500.440.440 
440  122-LL(2> 

T2-TU22  > 

445  IF ( LL ( 3 ) -XI ) 510.450.450 
450  LZ3-LH3) 

T3»T( LZ3 ) 

455  IF ( LL ( 4) -XI  >520.460.460 
460  LZ4-LU4) 

T4*T(LZ4) 

465  IFILL(S)-KI >530.470.470 
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470  LZ5«LLI5) 

T5«T(LZ5> 

475  IF(LL(6)-KI )540*480*480 
480  U6-LLI6) 

T6»T  <LZ6) 

GO  TO  600 
500  T2*1*E4 
GO  TO  445 
510  T3»1.E4 
GO  TO  455 
520  T4-1.E4 
GO  TO  465 
530  T5«1.E4 

GO  TO  475 
540  T6»1.E4 

600  01*AMIN1( T2*T3*T4»T5*T6) 

02«. 75*01 

KK«0 

TX»0. 

PT»0. 

TX2*0. 

TXPT-O. 

P2*0. 

DO  700  N«1 *500 
I F ( 02-T ( N ) ) 620 *620 *700 
620  I F ( Dl-T (N  >  >700*630*630 
630  TL ( N ) =ALOG( TIN)) 

TX-TX+TUN) 

PT  sPT+P ( N ) 

P2«P2+P ( N } **2 
TXPT»TXPT+TL(N)*P(N) 

KK«KK+1 
700  CONTINUE 
XtCK*KK 

B0« ( XKX*TXPT-TX*PT ) /<  XKK*P2~PT**2 ) 

AD«( TX-BD*PT)/XKK 
710  D3=EXP ( AD ) 

D6<=6  **D  1 

WR I TE ( 6*2220  >  01 
WRITEI6.2230)  03 

2220  FORMAT! 1H0.10X.33HAPPARENT  PUSiflVE  OURAT ION l MSEC ) ■ , 1F1Q.4// ) 

2230  FORMAT  I  10X.37HEXTRAPOLATE0  POSITIVE  01.RAT  ION  <M$EC  )  •  »1F10.4 ) 

RETURN 

END 

SIBFTC  SPRES 

SUBROUTINE  PRE >1 1  01 »T  »P  *MM) 

COMMON  X(4  »500)  »A(  10  >  *  PI  500  )  *UI  50)  *TU  500)  *PT  (  500 )  *  TPLOT  t  500 ) 
1PU500)  »0UMMV(50)  »IX<  2.500)  *T  I  500)  *U(500)  *TITLE<24)  ,01*0(500) 
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